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Effective classical Hamiltonians for spin systems: 
A closed form with quantum spin-wave effects 
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Thermodynamic properties of any quantum spin system can be described by the formally exact, 
although in general intractable, effective classical Hamilton function TL. Here we obtain an explicit 
form of TL which applies at T C JS 2 , where J is the exchange and S in the spin value, and 
incorporates quantum effects at the level of the spin- wave theory (SWT). For a quantum Hamiltonian 
H of Heisenberg form, TL is also Heisenberg but with a long-range effective exchange J*f , which is 
the price for including quantum effects. For three-dimensional magnets, classical SWT with TL yields 
the same results as quantum SWT for the original system, in the antiferromagnetic case with the 1/S 
correction to the ground-state energy. For nontrivial one- and two-dimensional systems, reduction 
of the problem from quantum to effective classical allows to apply such methods as classical Monte 
Carlo simulation and the 1/D expansion, where D is the number of spin components, etc. 

PACS numbers: 05.30.-d, 75.10.Hk, 75.10.Jm 



I. INTRODUCTION 

The partition function of any quantum spin system 
with a Hamiltonian H can be reduced to that of an ef- 
fective classical spin system with the Hamilton function 
Ti by rewriting it in the basis of spin coherent states 
|rij), thfyinit vectors rii serving as classical spin vectors 
in 7Y.BB'u The potential significance of the possibility to 
treat quantum systems as classical is quite apparent: In 
one and two dimensions one can perform classical Monte 
Carlo (MC) simulations, and use the classical-spin dia- 
gram techniqueaa or the 1/ZZ.expansion, where D is the 
number of spin components, BlU etc. However, since Ti 
incorporates all quantum effects of the system, in most 
cases its formally exact expression cannot be brought into 
a practically applicable form without approximations. 

In Ref. y a cumulant expansion (CE) was applied to 
semiclassical (S ^S> 1) one- and two-spin systems to ob- 
tain a tractable form of the effective classical Hamilton 
function TL. In Ref. || this approach has been generalized 
to many-body spin systems. The cumulant expansion 
yields TL as a series of non-Heisenberg multispin terms 
of increasing complexity, the zero-order term being the 
classical counterpart TL^ of the quantum Hamiltonian 
H. CE is an expansion in powers of JS/T with addi- 
tional terms of higher orders in 1/S. Thus it diverges 
for T < JS, where JS is of the order of the maximal 
spin- wave energy £j^ ax and is much smaller than another 
energy scale JS 2 which is of the order of the Curie tem- 
perature T c for 3d systems. Note that in the applica- 
bility region of the cumulant expansion the boson occu- 
pation numbers = l/(e^ Ek — 1) can be expanded in 
/3£k -C 1, that is, the system behaves classically with 
quantum corrections. Note that one should not confuse 



the cumulant expansion with the high-temperature series 
expansion (HTSE) which is an expansion in JS 2 /T and 
for S 3> 1 has a much narrower applicability range than 
the former. 

The presence of the two different energy scales in the 
problem suggests continuing the results of Ref. || into 
the region T J$ JS, which is the aim of this work. In- 
deed, at the convergence limit of the cumulant expan- 
sion, T ~ JS, there is already a strong short-range or- 
der (which is formed at T <C JS 2 ), i.e., fluctuations of 
spin vectors coupled by terms of the CE, relative to each 
other, are small. Thus the form of CE should simplify 
for T <C JS 2 , and one expects that the resulting series of 
effective Hamilton functions can be summed up to yield 
some TL that is valid down to T = 0. As we shall see, 
this works indeed for the ferromagnetic model, where the 
information contained in TL^ and TL^ of Ref. |^ is suf- 
ficient to guess all other terms of CE and to sum them 
up. 

For more complicated models, such as antiferromag- 
nets, higher-order terms of the simplified CE cannot be 
guessed from the first few ones, therefore one has to 
search for an alternative approach. The latter can be 
found. It turns out to be more satisfactory since one does 
not have to expand 7i in a power series and then sum it 
up again. Instead it makes use of the Holstein-Primakoff 
expansion of spin operators up to bilinear boson terms, 
which is familiar from the SWT. Thus our result for TL 
will be a good approximation in all cases where SWT 
works, i.e., for most of magnetic models except for those 
where quantum effects lead to nontrivial ground states. 

The paper is organized as follows. In Sec. [n] derivation 
of the general formula for the effective classical Hamilton 
function TL is given. Thereby the spin coherent states 
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are used and the results of the cumulant expansion of 
Ref. | are recapitulated. The CE for ferromagnets is 
simplified for T -C JS 2 and summed up to yield Ti which 
is valid down to T = 0. In Sec. Ill the Holstcin-Primakoff 



expansion is used to obtain the expressions for Ti for 
ferro- and antiferromagnetic models. In Sec. ^ possible 
applications of the results obtained are discussed. 



II. EFFECTIVE CLASSICAL HAMILTONIANS: 
FROM CUMULANT EXPANSION TO SPIN 
WAVES 

A. Main idea and basic relations 

Thermodynamic properties of any quantum spin sys- 
tem with a Hamiltonian H can be described in an alter- 
native way via the effective classical Hamilton function 
Ti which can be constructed using spin coherent states 
|n,*) on lattice sites z.Eratl The vectors rij parametrizing 
these maximum-projection states serve as classical spin 
vectors in Ti, and the partition function Z of the quan- 
tum spin system can be calculated as an integral over 
drij. The derivation of this effective classical formalism 
for spin systems is based on the representation of the 
unit operator in the (overcomplete) basis of spin coher- 
ent states 



25+1 
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dn |n) (n| 



(1) 



For a single-spin system, the trace of an operator A over 
any complete orthonormal basis \m) becomesH 



tr A — ^2(m\A\m) = ^ + 1 / dn V7m|A|n)(n|', 



25 + 1 

47T 



dn (n\A\n) 



(2) 



Similarly, the partition function for a many-spin quantum 
Hamiltonian H can be written as 

Z = (™+±Y I ^^({n^le-^Kn,}), (3) 
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i=l 



where (3 = 1/T. It has the same form as the one for clas- 
sical systems, provided one defines the effective classical 
Hamilton function Ti by the formula 



e -/m( { n i}) = ({n . } | e -M| {n . }) . 



(4) 



The function Ti evidently depends on temperature, 
thus calculation of physical quantities should be done 
with care. For instance, the internal energy is given by 



U 



dlnZ _ / d(pH) \ 
d(3 ~ \ ap / 



(5) 



where (. . .) denotes a classical thermal average with Ti. 
One can see that Ti* = d(/3Ti)/d(3 ^ Ti. The same care 



should be taken in calculation of the magnetization and 
correlation functions (see Ref. ||) . 

The matrix element over the spin coherent states which 
defines Ti in Eq. (|4|) can be calculated explicitly-only 
for the simplest models such as one spin in a fieldfl For 
many-body systems such as the Heisenberg model 



H 



(6) 



one has to resort to approximate methods. A plausible 
idea is to solve the problem perturbatively in 1/5 in the 
quasiclassical case 5> 1. In the classical limit, 5 — > oo, 
one obtains the classical Hamilton function 



i ij 

as was rigorously proved by Lieb.S Here 



5Hj, 



5 Ji 



(7) 



(8) 



are the reduced magnetic field and the exchange interac- 
tion, respectively. 

Quantum corrections to Eq. (0) were calculated in Ref. 
|§| with the help of a cumulant expansion which yields a 
series of multi-site non- Heisenberg Hamilton functions of 
increasing complexity. These expressions are too lengthy 
for reproducing them here. The small parameter required 
for the validity of the cumulant expansion, in addition to 
the obvious 1/5 -C 1, is 



JS 1 



(9) 



where J is the coupling constant when only nearest neigh- 
bor interactions are taken into account. The inequality 
breaks down at a temperature T ~ JS which is of the 
order of the maximal energy of spin waves. On the other 
hand, a strong short-range order sets for T < JS 2 , which 
is much higher than JS for S>1. Thus for 5 ^ 1 the 
results of the cumulant expansion should simplify in the 
temperature range JS «T« JS 2 , and one expects that 
CE can be summed up in this range to yield an expres- 
sion for Ti which is valid down to T = 0. This is the idea 
of the present work which in the following subsection is 
presented for ferromagnets. 



B. Effective low-temperature Hamilton function for 
ferromagnets 

The non-Heisenberg terms of the cumulant expansion 
for the Hamilton function Ti (see Ref. ||) can be simpli- 
fied for T <C J = JS 2 since there is strong short-range 
order in the system in this temperature range. For the 
ferromagnetic model which we will deal with in this sub- 
section, the scalar product rii • rij is close to unity. Thus 
one can neglect the terms quadratic in 1 — rij • rij and 
simplify 
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(rij • nj)(nj ■ n ; ) = [1 + (n; ■ n,) - 1][1 + (n, • n ; ) - 1] 
-1 + (n,- • jLj) + (nj • n ( ) (10) 

With the help of this formula, the first-order Hamilton 
function in zero field, 

H(1) = E 4^K^ • n,) - (n 2 • nj )(n 3 • n,)] (11) 



which was derived before [see Eq. (20) of Ref . || , simpli- 
fies to 



with 



V iP = - 2 + -~- JilJlj- 

Jo , 



(12) 



(13) 



The second-order terms of the Hamilton function of Ref. 
H can be processed in the same way. The resulting ex- 
pression for TL up to the second order reads 



with 



and 



,12. , _j_ /^0y(l) _ (^O) 2 T/ (2) 



$ ) V = J, 



tj "y 1 25* *J 6S* 2 J ' 



(14) 



(15) 



3 1 

V%j = ^ij^O — 3Jy + -=- ^ J; Jy =r ^ JilJlnX 

Jo , 



J In 



(16) 



In Fourier representation Eq. ( |l5| ) takes the form 



$ k 2) = Jo - (Jo - Jk) + t^(Jo - J k ) 2 - ^(Jo A ) 



2S ,( 



65 2 



(17) 



It is not difficult to guess the form of TL containing all 
orders of PJo/S in the cumulant expansion at T <C J — 
JS 2 . The result reads 



with 



$ k = Jo - TS[1 - exp(-/fe k )], 



(18) 



(19) 



where e k = (Jo — Jk)/S — S(Jq — J k ) is the spin- wave 
spectrum of the ferromagnetic model. 

One can generalize Eqs. ([l8]) and ( |l9| ) for the ferromag- 
netic model by adding the effect of a homogeneous field. 
This is done by applying the same method to the field 
terms of the cumulant expansion of Ref. ||. The result is 



W = Ho-h off -^n J -i^J^n J .n J , (20) 



where 




with the effective field and the exchange interaction given 
by 



h cff _ TS-(1 



-PH\ 



jf = TSe~ l3ek . (22) 

The spin-wave spectrum is here of the form 

£ k = (h + Jo - J k )/S = H + S{J - J k ). (23) 

One can check that in the limit S — ► oo the effective 
Hamilton function of Eq. (|2^) reduces to the classical 
form of Eq. ^ . The minimum of Eq. (^p| ) is attained for 
all spin vectors directed along h and is equal to 



Eq = -Nh - 



NJ 



(24) 



i.e., it coincides with the ground-state energy of the clas- 
sical model. 

It should be noted that although the form of TL given 
by Eq. (|o|) has been derived for T <C JS 2 , it yields qual- 
itatively correct results for the thermodynamic functions 
in the whole range of temperatures. For T 3> JS 2 , devia- 
tions from the exact high-temperature asymptotes are or 
the relative order l/S, whereas the leading classical term 
is recovered (see Appendices A and B for toy models). 

Analogous calculations can be done for the antiferro- 
magnetic model on bipartite lattices etc. For the bipar- 
tite AFM model in zero field one obtains Eq. (113) with 



(2) --Jk + ^(Jo + Jk) 2 



6S 2 



(Jo 



J k ) 2 (J + 2J k ). 

(25) 



Unlike the case of Eq. (|l7j), here it is difficult to guess 
a function which could yield an expansion of this form. 
Thus another approach to this problem is needed. The 
possibility of constructing the low-temperature form of 
the effective classical Hamilton function for the ferromag- 
netic model is an indication that an analytical method 
must exist which is not based on summing up the cu- 
mulant expansion and which works for both ferro- and 
antiferromagnets . 



III. CLOSED-FORM QUASICLASSICAL 
HAMILTONIANS WITH SPIN WAVES 

A. A boson formalism 

Returning to the basic formula for the effective Hamil- 
ton function, Eq. (M), we note that in many works the 
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introduction of spin coherent states for large-spin sys- 
tems is followed by the assertment that such spins behave 
quasiclassically. In that case one can simply replace the 
spin operator by Siii, which obviously leads to the 
classical Hamilton function. Such an approach, however, 
misses quantum fluctuations of spins which are respon- 
sible for Si 7^ Sn.i and therefore possible corrections to 
the classical results. To take quantum effects into ac- 
count, we will proceed as follows. First, we express spin 
operators in the coordinate system defined by the vector 



where 



Si — Si Z Tii + Si+n 



Si-XL 



i + 



n i± = (n„ ± in iy )/2. 



(26) 



(27) 



The vectors \\i X and rij a are unit vectors which form, 
together with = n^, an orthogonal coordinate system. 
Then we use the Holstein-Primakoff boson representation 

S.f = S - a+a l: S+ = (2S - afa,) 1 / 2 ^, 

Sr = a+(2S-a+a l ) 1 ^. (28) 

Here we expand square roots and keep in H only up 
to quadratic terms in boson operators. A priori, this 
procedure is justified for S > 1. In fact, however, it is 
similar to the usual approach made in the linear spin- 
wave theory, and in most cases it does not require large 
S. The latter thus becomes 



(29) 



Here TL^ is the classical Hamilton function given by Eq. 
(0), 



tf (1) = -A/~£[(h 



MF 



n, 



where 



,MF 



(30) 



(31) 



thus the scalar products • rr, are close to their zero- 
temperature values n j ■ ir, . For this reason, we will re- 
place in Eq. (^) rij + • n,-± by rij+ • nj-±, etc. The advan- 
tage of it is that the n i+ ■ Hj± arc translationally invari- 
ant and this allows us to rewrite the bosonic quadratic 
form in Fourier space as a sum over a single wave vector. 
This greatly simplifies the solution. We will also con- 
sider two toy models where the problem can be solved in 
a general, i.e., non-simplified form. These are the spin- 
in-a-field model (Appendix A) and the two-spin model 
(Appendix B). 

The choice of the transverse vectors and in Eq. 
( p7| ) is not unique. One can rotate them in the plane 
perpendicular to n^, which does not change the final re- 
sults for physical quantities. The quantities n i+ • TiJ± in 
Eq. ( |32| ) can be specified already at the present stage by 
an appropriate choice of and n^. In most of mag- 
netic models at zero temperature the vectors lie in a 
plane. In that case it is convenient to choose equal 
for all sites i and perpendicular to this plane. It follows 
automatically that = [n^ x n^]. This yields 



1 



and thus 



n l± • n jT = (1 + m ■ nj)/4 
n,± • n J± = (1 - n~TT~)/4. 



(33) 



(34) 



Evidently the definitions above remain valid if the spin 
vectors in the ground state are collinear, as for ferromag- 
nets or antiferromagnets in zero field or in strong fields. 

Using the formulas obtained above in Eq. ( |32j ) and 
Fourier transforming one obtains 



H 



(i) _ 



- X(.9kak + g^) 

k 



with 



3k 



(35) 



(36) 



is the classical molecular field, and 



h (2) = -£ X ^[( ni + ■ n o-) a i a 3 + ( n i+ ■ n J+y 

+ h.c.] + i^(hf F . ni )a+ a4 . 



(32) 



The problem of finding the effective classical Hamilton 
function Ti for a spin system, Eq. (|]) , has been reduced, 
in this approximation, to calculating the average of an 
exponential of a quadratic form of boson operators over 
the bosonj-vacuum. The latter can be done analytically 
in generalEl but the result for a many-spin system is so 
complicated that it has no practical significance. On the 
other hand, in the most interesting temperature range 
T -C J = JS 2 there is a strong short-range order, and 



and 

k 

where 



and 



A k aja k + 7; B k( a k a -k + n - c - 



■J i 



13 j 



(37) 



A k = SJ' - S( J k + J k )/2 + H • n, 
B k = S( J k - J k )/2, 



(38) 



jL = j2 e iMn-r j ) J .. W rn- j 



(39) 
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In Eq. (^) it has been assumed that in the ground state 
H ■ rii does not depend on i. This is true for ferromagnets 
and bipartite antiferromagnets. 

Now we are prepared to perform the main step of our 
method and to calculate the matrix element over the bo- 
son vacuum: 

H = H {a) - T\n(O a \ exp[-0(6W + # (2) )]|0 a ), (40) 

cf. Eq. (Q). Since bose operators with different k com- 
mute, [ok, Oq] = A(k— q), and -ff' 2 ' couples bosons with 
k and — k, the exponential in Eq. (40) can be factorized 
in a product of exponentials containing different pairs of 
ak and a_ k . Thus our many-body problem reduced to a 
problem of calculating the vacuum matrix element of the 
exponential of a bilinear form in bose operators a = a k 
and b = a_ k (see Appendix C). The result has the form 

H = + E Q - ^ + <?_k)G?k + g*_ k )F+ 

k 

-\T,^-9-^-9- k )F k (41) 



where e fe = a/<5J" 

F± = 1 
k A u ± S k 

and the quantity 



B\ is the spin-wave spectrum, 
2T 



A k ±B k + e k coth(/3e k /2) 



(42) 



Air 

Tin ( cosh(/3e k ) + — sinhf/fek 



(43) 

will be shown to be related to the spin-wave quantum 
correction of the ground-state energy for antiferromag- 
nets. 

Using Eqs. (|36|) and (|27|) one can rewrite 7i as 



n = n 



E c 



- G y F~] 



(44) 



where 



Gf, = (hT-n 



and 



{hi ■ ii„,m1i. 



)(hf F 

MF 



n jy ) 



(45) 



(46) 



The molecular field h^ F given by Eq. ( pl| ) is a linear 
function of {n;}. Therefore the quasiclassical Hamilton 
function Ti. of Eq. ( [flf ) contains four-spin interactions 
and is thus of a non-Heisenberg type. Eq. ([fZJ) can be 
simplified to a Heisenberg-type Hamilton function in the 
most interesting temperature range T -C J = JS 2 by 
using the fact of a long- or strong short-range order at 
these temperatures. 



Before proceeding to this, let us consider on passing the 
quantum corrections to the classical Hamilton function, 
W°\ at high temperatures. To the leading order in 1/T 
one obtains for Eq 



16S 2 



(47) 



Using i 7 ^ = /3/2 to leading order in as well as re- 
lations from the Appendix of Ref. ||, one obtains for the 
remaining quantum corrections in Eq. (Ell) 



« (1) = -§5>k 5k = -|;£[h: 



(48) 



[The same result follows if one keeps only the linear boson 
terms in H, i.e., by setting A k = £? k = in Eq. (pil|)]. It 
can be seen that coincides with Eq. (20) of Ref. || 
without the last small term of order (3/S 2 . 

To obtain a quasiclassical Hamilton function which 
is suitable in the low-temperature range, one has to 
make our approach self-consistent by applying a small- 
fluctuation approximation, as was already done above in 
order to simplify Eq. (^) . For three-dimensional systems 
with long-range order at low temperatures, this would be 
a straightforward task since one can write n; = + 6n.j 

and expand in small Siu around the long-ranged n,- '' . For 
one- and two-dimensional systems without long-range or- 
der, there is no nf\ yet deviations from the short-range 
order are small. Exploiting this smallness requires a more 
sophisticated and general approach which we are going 
to present now. 

First, in Eq. (E5L) it is convenient to replace 



,MF 



h™ F - h 



MF,0 



n. 



(49) 



where h i ' is the zero-temperature value of /i^ F . Since 
in the ground state spin vectors are pointing along the 
molecular field, Qi deviates from zero by virtue of fluctu- 
ations only. Thus in Eq. ( |45| ) which contains Q quadrati- 
cally, one can neglect fluctuations of vectors etc. Sec- 
ond, we expand Q; in the local basis of site i as 



(Qi • n,-)nj + (Qj (50) 



and perform the scalar product Q; • Qj, then we simplify 
it with the help of Eq. (|33|). Using the resultin g fo rmula, 
we eliminate (Qi • Ui X )(Qj ■ nj x ) in Gfj in Eq. (J45|), after 
which Eq. (44) becomes 



U = W + H" 



where 



and 



(51) 



(52) 
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H" = 



^+[n; • njiQi ■ ni)(Qj • n,) 
•ni-n iv (Qi • n i )(Q j • n jy ) 
n i yn j (Q l ■ ni„)(Qj • n 3 -)]}. 



(53) 



The advantage of representing 7i in the above form is 
that H" vanishes for models with collinear ground state. 
Indeed, in this case • n.j y = 0. Further, the quantity 
Qi • n.j becomes quadratic in spin deviations and thus the 
term with (Q^ • rii)(Qj ■ rij) is of fourth order and can be 
neglected. The remaining term in TL" vanishes because 



of the factor (F- 



rii • njF^) which can be shown to 



be zero for collinear states. For ferromagnets, one has 
rii ■ Hj = 1 and = F^j , while for antiferromagnets in 
zero field one has n; • n , = ±1 and F^ = ii 7 !" depend- 
ing on whether i and j belong to the same or different 
sublattices, respectively. Thus Ti" is nonzero only for 
the bipartite antiferromagnetic model in the field range 
< H < Hsf, where Hsf = 25 Jo is the spin- flip field, 
as well as for more complicated models which are not the 
subject of this paper. Next we will consider the ferromag- 
netic and bipartite antiferromagnetic models separately. 



B. Quasiclassical Hamilton function for ferromagnets 

For the ground state of the ferromagnetic model one 
has (h • nj) = h and (n; • n,) = 1. Thus the parameters 
defined in Eqs. ( |38|) and @) simplify to 



J k — Jk, i?k — 

A* = S(J - Jk)+# = £k, 

which results in Eq — and 

= = = ^{e- 0ek + /fek - 1) 
Thus Eq. (|5l]) reduces to 



where 



Q 4 = hf F -(h + J ) ni 



(54) 



(55) 



(56) 



(57) 



Recalling the definition of h^ F given in Eq. (|3l"|), one can 
write 



-ik(r« — Yj ) 



k ijlm 

x [6u(h - /in,) + (J u - S a J )n z ] 

• [<5jm(h - hnj) + (Jj m ~ 6jmJo)n m ] 

W (0)_ JLy Fk y e-iHn-rj) 
2N ^ ^ 

k ij 

x [H - £ k ni] • [H - e k nj] 



W (o)_^! Fo + i/Fo ^ (h . ni) 



S 



(58) 



where [e^jF]k]y denotes Fourier transform of sf.F^. With 
the help of Eqs. (Q) and ( |55|) this result can be brought 
into the final form of Eq. fl2C|). 



C. Quasiclassical Hamilton function for 
antiferromagnets 

For the antiferromagnetic model it is convenient to re- 
place J — J in all the formulas above. The canting 
angle of a classical bipartite AF model in an applied field 
is given by 



ana 



H 

2SJ 



h 
2X 



(59) 



for ft, < 2 Jo and x = 1 for h > 2 Jo. Thus in the zero- 
temperature equilibrium state one has 



H rij = xH 

xtTnj = x 2 + (I ~ x 2 )e- lb{r 



(60) 



where b is the inverse lattice vector (for the simple cubic 
lattice b = {7r,7r,7r}). From Eq. tt 
(2x 2 — 1) Jk, so that 



one obtains J£ = 



A k = S{J + x 2 J k ) 
B k = S{1 - x 2 ) J k 
A k + B k = S{J + J k ) 

A k - B k = S{J - J k + 2x 2 J k ). (61) 
The antiferromagnetic spectrum has the form 

e k = Sy/(J + Jk)(Jo- J k + 2a; 2 Jk). (62) 

In fact, this expression describes the two spin- wave 
branches of an antiferromagnct in field, i.e., the gapless 
£ik = £b-k = 5a/(J - Jk)(Jo + Jk - 2x 2 J k ) and the 
gapped £2k = £k one. Whereas £2k describes spin mo- 
tions within the plane specified by the equilibrium ori- 
entation of the two sublattices (this motion changes the 
Zeeman energy), the gapless £ik accounts for the motions 
out of plane, which cost no Zeeman energy. The quantity 
Qi defined by Eq. (W9) is given by 



Qi = hf IF - J n, 



(63) 



[cf. Eq. Q57|)]. Now one can process W of Eq. (|2J) in 
a way similar to that of ferromagnets in the preceding 
subsection, which yields 



W = H (0) 



E Q - ^^F+ + 2SJ F+ £(h • n, 



IS 



[(A k +5 k ) 2 F+]«K-ni) 



(64) 
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[cf. Eq. dH]. With the help of Eqs. @) and @ this 
equation can be brought into the final form 

H> = Wo + E Q h cff • + ~ E -C'»< • n ^ (65) 

i " ij 

where Eq is given by Eq. (^) , 
7Yo 



where has vanished. Now one can compute Qi • n.^ 
with the help of Eq. (f70f) and the formula 



x = - y^fe Jo 
i 



Ju]{ni 



V). 



(74) 



N h(h-h° B ) | Jp 
4J 2 



The result has the form 

ihvi 

ly ~ 



Q 



- — - 



X 



s/T 



x- 



(66) 



i 



The effective field and the exchange interaction have the 
form 



and 



,cff 



JcB,± 



2STH 



(75) 

The expressions obtained above can be substituted into 
Eq. (pis]), which after some algebra yields H" in the 
anisotropic Heisenberg form 



A + B a + e coth(/3e /2) 

2ST{A k ± gk) 
v4 k ±Bk + £k coth(/3e k /2) ' 



(67) 



(68) 



H" 



-^Gyflli-X-I^Ilj-X 



where 



(l-^)jf- + + x a ^+-J£_ k 
x 2 (l — x 2 ) 



TCff, 



(76) 



(77) 



Expanding jjf '+ for /3e k ~ SJ/T <C 1 one recovers Eq. and jcff,± are defined b E ^ ft can bc shown that 



(||) in the zero-field case. Clearly, it is very difficult to 
guess that the series in Eq. ( p5| ) can be summed up to give 
Eq. ([38]) . It can be checked that in the classical limit Eq. 
( pq ) goes over to the classical Hamilton function 1iS \ 

We calculate next U" of Eq. (||). To this end, we 
define 



2Jr 



(69) 



[cf. Eq. (|59|)] and use the explicit definition of the trans- 
verse vectors 

_ [m x x] Jhr . _ [m x x] JhT . 



Gk remains finite in both limits x — ► and x — > 1, and 
therefore, 7i" vanishes in these limits. 

Let us consider the minimal value Eq of the antifer- 
romagnetic quasiclassical Hamilton function, Eq. ([5]]). 
Minimization of Eq. (|5l]) yields the value of the canting 
angle which is the same as for the underlying classical 
model [cf. Eq. ©]: 

h cS h 
x = -^ w -r = -^. (78) 

For this x the term 7i" vanishes. Using Eq. (^o|) and 
noticing that j£ :+ = one obtains 



xyl — x 2 



n ly = [m x n ix ] = 



TO-iX — X 

X\/l — X 1 



En — E, 



(70) 



ci 



E c 



Eci — — 



Nh 2 NJ 
4Jo~~^~' 



(79) 



(remember that we neglect fluctuations in rij as all fluc- 
tuations are accounted for by Qj). With these definitions 
and Eq. (|6C|) one has 



Here Eqi is simply the ground-state energy of the clas- 
sical model, as in the ferromagnetic case, cf. Eq. J24|). 
The quantity Eq given by Eq. (^3|) is of quantum origin 
and can be rewritten in the form 



Vl-i 2 [l-e' b(ri ^]e' 



br ; 



E c 



AE -T^2 ln(n k + 1) + | ^ ln 



• n iy — —xyl — x 



(71) 



(ST) 2 



feff,+ TeS,- 
J k J k 



Next we simplify the quantity • which can be rewrit- 
ten in the following way 



(80) 

where n k = 1/ (e^ £k — 1) is the boson occupation number 
and 



Q 4 • Jon, = -Q» • (Qi + 2J m) - ^Q 2 



1 



A£b = -±$>l k 

k 



£k). 



(81) 



2(hf lF -J 2 ). 



(72) 



Here we neglected the quadratic term Q 2 . Substituting 
for h^ F the expression given by Eq. (|3l|) one obtains 



Q 4 • rij = -2 E -M n i ' x - x2 )i 



(73) 



is the first-order correction to the quantum ground-state 
energy. In the limit T — > one has £7o = Sci + A£o. 

Again, we stress that, although derived for T <C J,!? 2 , 
our effective quasiclassical Hamilton function Tl leads 
to qualitatively correct results in the whole temperature 
range, reproducing the leading classical terms in the high- 
temperature asymptotes of physical quantities. 
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IV. ILLUSTRATION: SPIN- WAVE THEORY 



V. DISCUSSION 



For an illustration of the formalism developed above 
we consider a three-dimensional magnetic systems for 
which we know the thermodynamic functions derived by 
the standard spin-wave theory. We want to demonstrate 
that classical SWT for the effective Hamilton function TL 
yields quantum results. Let us do it for the antiferromag- 
netic model. Starting from the zero-temperature ordered 
state one can write 

(0) i X 

Sn t - ngxi + nfJ Vl - (z? + yf )/2, (82) 

where Xi, yi <C 1. Then TL' of Eq. ( |65| ) becomes a bilinear 
form of spin deviations Xi, yt 

H = E Q + - 2^ h i ' (Xi+Vi) 

i 

+ Jif' + { x i x i + n ly • ri jy y iyj ) (83) 

where ■ n,-„" = • n j . Furthermore the energy Eq is 
given by Eq. (I79j), and the effective molecular field is 

h MF,eff,0 = h eff_^ i ^,+ n (0) - (g4) 
3 

For an antiferromagnet in an external field h < 2Jq one 
has hf F '° ff '° = 0, as follows from Eq. @ and Jf' + = 0. 
For TL" of Eq. @ we obtain 

H " ~ ~\ H G' ij (x-n iy )(x-n jy )y 4 ^. (85) 

Adding Eqs. (f|) and © and using Eqs. © and (0) 
yields TL in the harmonic approximation 

ij 

After calculating the partition function and using Eqs. 
( f79| ) and ( ]80| ) one obtains the formula 

In Z = const - /3(Ea + AE Q ) + ^ ln(n k + 1) (87) 

k 

which is nothing else but the result of the spin- wave the- 
ory for quantum systems. 

For the ferromagnetic model the considerations are 
much simpler. One has h^ F:C '° = TS and the harmonic 
form for Tt reads 

TS 

H = E + —Y,[^-- e'^ijixiXj + y iVj ). (88) 

It yields an expression for In Z similar to that of Eq. ( p^ ) 
with Eq and £k appropriate for ferromagnets. 



We have obtained an effective classical Hamilton func- 
tion H. for ferromagnets and bipartite antiferromagnets, 
which can be used to compute thermodynamic quantities 
of the original quantum spin system with the Hamilto- 
nian H. Quantum effects are taken into account in TL 
at the level of the linear spin-wave theory, which means 
that the results can be applied in all cases where the 
quantum corrections to the ground state are small and 
can be treated perturbatively. Although the results for 
TL have been obtained for T <C JS 2 , the thermodynamic 
quantities computed with TL remain qualitatively correct 
in the whole range of temperatures. Indeed, since TL re- 
covers the classical Hamiltonian in the limit S — * oo, the 
main classical terms in the high-temperature asymptotes 
are correctly reproduced and the errors are of relative 
order 1/S. 

Our results can be applied to one- and two-dimensional 
spin models which are characterized by the absence of 
long-range order. Here a reduction of the quantum model 
to a classical one is an important step in simplifying the 
problem. However, the main difficulty in low dimen- 
sions is due to long wavelength fluctuations which pre- 
vent ordering at any finite temperature, and this feature 
is characteristic to both quantum and classical systems. 
Thus solutions of effective classical models for quantum 
systems are not straightforward. They require applica- 
tion of numerical methods such as Monte Carlo simula- 
tions or apejepximate analytical methods such as the X/D 
expansion,Lli3 where D is the number of spin components. 

The latter is very efficient for systems such as kagome 
and pyrochlore antiferromagnets which have a degener- 
ate classical ground state and remain a classical spin liq- 
uids (i.e., they do not show any extended short-range 
order) down to zero temperature. For these two models, 
even a D — > oo approximation yields excellent results for 
thermodynamic quantities, as compared to the numerical 
data. Studying quantum effects in these systems with the 
help of effective classical Hamilton functions in combina- 
tion with the 1/D expansion seems to be very promising. 
However, for these models the form of TL will be different 
and it is still to be determined. 

In concluding, one should mention an alternative ap- 
proach to the derivation of effective classical Hamilton 
functions which was developed by the Florence group 
(see Refs. |^, and references therein), which, in combi- 
nation with MC simulations, proved to be very efficient 
for two-dimensional antiferromagnets. The main differ- 
ence between our and the Florence forms of TL is that 
our approach leads to a long-range effective exchange in- 
teraction Jfj with a diverging range of interaction at 
T — » [see, e.g., Eq. (]22j)], whereas the expression for 
Jfj of the Florence group remains restricted to nearest- 
neighbours for the quantum models with nn interactions 
while only the magnitude of J is renormalized by quan- 
tum effects. The latter makes it possible to solve one- 
dimensional ferro- and antiferromagnetic models by using 
the known solution for classical systems. This leads to 
results that are in a good accord with the numerical data 
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for the appropriate quantum models In distinction, an 
analytical solution for the classical magnetic chain is not 
available for long-range interactions. 

While the long-range form of J?P resulting from our 
method might be considered as a disadvantage from the 
practical point of view, we should stress that it ap- 
pears in a natural way and has a physical explanation. 
The integrals over the Brillouin zone in the usual SWT, 
such as the magnetization correction in ferromagnets 
Am = — Jc?knk, are at low temperatures restricted to 
the small-k region satisfying /3ek < 1 because spin waves 
with larger k are frozen out. In our effective classical 
formalism, the same cut-off is ensured by the effective 
interaction, i.e., by J£ ff = TSe~^ £k for ferromagnets. 
Clearly localization of J£ ff in the small-k region at low 
temperatures means that the interaction is long ranged 
in the coordinate space. 
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In the limit T ^> H it reproduces the first three terms of 



the cumulant expansio 



H = -S(H-n)^[Hxnf 

+ ^(H.n)[Hxn] 2 + 0(/3 3 ). 



(A7) 



Unfortunately, at low temperatures Eq. ( |A6| ) leads to di- 
vergencies in the thermodynamic properties because the 
spin tends to direct itself against the field and to cre- 
ate a large number of bosons. This takes place since we 
have kept only quadratic boson terms for spin operators. 
To avoid such a divergence, we linearize H. with respect 
to fluctuations near the zero-temperature state with the 
spin directed along the field. The linearized Hamiltonian 
has the form 



H = -SH + TS(1 



1 - 



H n 



H 



(A8) 



Now one can easily calculate the partition function 

' 1 _ e -2S(i-e-f> H y 



In Z = const + j3SH + In 



1 - e-f> H 



, (A9) 



APPENDIX A: ONE SPIN IN A FIELD 

In this Appendix we shall apply our approach to the 
most simple model, i.e., one spin in a field: 



H 



H S. 



(Al) 



In boson representation in the basis of spin coherent 
states the spin operator takes the form 



(A2) 



S = (S - a + a)n + V2S[n + a + + n_a] 



[cf. Eqs. (|26|) and @]. The Hamiltonian of Eq. ([At]) 
thus becomes: 

H = —(S — a+a)(H • n) - V2S[(H • n+)a+ + h.c.]. 

(A3) 

For this model we are able to calculate the effective 
classical Hamilton function 



H = -Thi(O a \exp[-0H]\O a ) 



(A4) 



at all temperatures, i.e., without assuming small fluctua- 
tions of the vector n. The corresponding matrix element 
over the boson vacuum is calculated easily resulting in: 



(0 o | exp[-/3(Aa + a - g*a + 
= exp 



ga)]\o a ) 



L9 {e -0A + /3A _ 1] 



(A5) 



Therefore, the effective Hamilton function becomes: 
ST [H x n] 2 



W = -5(H-n)- 



2 (H-n) 



[ e -m-n + ■ n - 1] 



(A6) 



which is well behaved at low temperatures. It correctly 
reproduces the exponential behaviour of the specific heat 
C . For large values of spin the deviation from the exact 
result is exponentially small: 



^exact — c 



25 



e 2S - 1 



(A10) 



This approach gives also correctly the leading classical 
term of the specific heat at high temperature 



C 



S(S + 3)H 2 



c _, a 32±l>^ (A11) 



^rp2 5 exact 2>T 2 

but there is an error in the 1/S quantum correction. 

APPENDIX B: TWO SPINS 

Let us consider the model of two ferromagnetically cou- 
pled spins S: 



H 



-JSi -S 2 . 



(Bl) 



First we rewrite the spin operators in terms of coherent 
states n, m and two bosons af, 

Si = (S - a^ai)n + V25[n + a+ + n_Oi] 



S 2 = (S - a+a 2 )m + V2S[m + aJ + m_a 2 ], (B2) 
after which the Hamiltonian takes on the form: 
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H = H (a) +H^+H^ 

H {1) = -jV2S^[(n+ • m)a+ + (n ■ m+)a+ + h.c] 
H {2) = JS(n ■ m)[o^Oi + a+a 2 ] 

— J5[(n+ ■ m_)a^a2 + (n+ ■ m + )afa 2 + h.c. 



(B3) 



After a canonical transformation to new boson opera- 
tors 



a+ = (b+ + b+)/V2, 

the Hamiltonian becomes 

H = H {0) + Hi+H 2 
R 

Hi 



(bt-bt)/V2 (B4) 



#2 

with 



-(#*&+ +gb 2 ) + 



Ab+b 2 



(i4 



A = J Ai 

2 y 



3n • m), 

g = US 3/2 \[n x m 



B = J Ai 

2 y 



bl) 



n • m 



(B5) 



(B6) 



Since H± and H 2 commute one can perform indepen- 
dent vacuum averages for the two kinds of bosons. For 
Hi one can use the formula: 



(0 bl \e x ^ +b ^\0 bl ) = 



(B7) 



whereas the average with H 2 can be calculated by a 
method similar to that presented in Appendix C. This 
leads to the effective Hamilton function of the form: 



H = W<°> 



B 
~2 



T 



ln(l -(3B)+E Q -g*gF-, (B8) 



where F is defined by Eq. (M2|) without the index k. The 
function Eq is given by Eq. ( [f3[ ) without summation over 
k, and £ = \J A 2 — B 2 . 

Although this Hamilton function gives correctly the 
first two terms of a cumulant expansion, it breaks down 
at low temperatures. Like a single spin in a field, at low 
temperatures the two spins tend to orient antiferromag- 
netically with respect to each other and create a large 
number of bosons, which results in a divergence of the 
partition function. In order to avoid this divergence, one 
has to expand TL up to the quadratic terms in the relative 
deviation of the vectors n and m, that is, one has to keep 
1 n • m and to neglect (1 — n • m) 2 etc. As a result the 
effective Hamilton function becomes: 



i _ -2/3JS 

H = -JS 2 + TS (1 n • m) 



which leads to the partition function 



InZ = const 



In 



1 _ e -2/3JS 



(B9) 



(BIO) 



At low temperatures one obtains an exponentially small 
specific heat, which differs from the exact result by an 
amount of order 1/ S: 



C 



AJ 2 S 2 



,-2(3JS 



J^2 



s 



U 2 S 2 
2^2 ' 



-2/3JS 



e s -l 
1 



45 + 1 



(Bll) 



At high temperatures our approach gives correctly the 
leading classical term: 

J 2 S*(S + 6) „ J 2 S 2 (S + l) 2 

U — r^rp2 ' ^CXact — ^rp2 ' l 131 ^ 

whereas there is an error in the 1/5* correction. 



APPENDIX C: VACUUM AVERAGE OF THE 
BILINEAR-BOSON-FORM EXPONENTIAL 

Here we calculate the exponential of a bilinear form in 
boson operators over the boson vacuum: 

f{x) = (0|e-^|0), (CI) 
where the operator H has a form: 

H = A(a + a + b+b) + B(a + b + + ba) 

+ 3> + + 9* b b + + g a a + g b b. (C2) 
First let us consider the following operators: 

a+{x) =e- x "a+e x " 

b{x) = e- xA be xA . (C3) 

Differentiating with respect to x, we find that a + (x), b(x) 
satisfy the system of equations: 



da + {x) 

dx 
db{x) _ 
dx 



-[H,a + (x)} = -Aa + (x) ~ Bb(x) - g a 
- [H, b(x)} = Ab{x) + Ba+ (x) + g* b . 



Integrating it, one obtains 

a + (x) = ua + + vb + Zi 
b(x) = u'b + v'a + + z[* , 

where zi = wig a + w 2 g* b , z[ = w[g a + w' 2 g^, 



( U V ^ 






-A 






= exp 


[.( 


B 





and 



-wi w 2 



u — 1 v 
v 1 u' - 1 



-A -B 
B A 



(C4) 

(C5) 
(C6) 

(C7) 



The explicit forms of u, v, wi and w 2 are 

, . . . sinh(a;e) 
u = cosh(xe) — A 

£ 

sinh(x£) 
v = —B 



Wl 



w 2 



cosh(a;£) — 1 sinh(x£) 



-B 



cosh(a;£) — 1 



(C8) 
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with e = V 'A 2 — B 2 . The quantities v! , v' , w[ , w' 2 are 
deriv ed f rom u,v,wi,w 2 by changing the sign of x in 
Eqs. (Cg). For the pair of operators b + (x) and a(x) one 
can obtain the similar to Eqs. (C5) formulae: 

b + (x) = ub + + va + z 2 
a(x) = u'a + v'b + + z' 2 * (C9) 

with z 2 = wiffb + u>25*, 4 = + w 29* a - 

Before calculating f(x) it is useful to introduce an aux- 
iliary functions f a {x) and fb{x) 



f a (x) = (0|e-^a+|0) 
f b (x) = (0|e-**6+|0). 



(CIO) 



After carrying the ope rators a + and 6 + past e xH in the 
last terms using Eqs. ( |C5| , |C9|) one finds: 



f a (x) = (0|a+(a:)e-* ff |0) = «/ 6 *(a;) + zi/(z) 
f b (x) = (0\b+(x)e- xfI \0) = + z 2 f(x) (Cll) 

and, therefore, one can express f a (x) and fb(x) through 



fa(x) 

fb(x) 



1 — ir 

Z2+vzt m 



(C12) 



1-v 2 

In order to find f(x) we perform the following manip- 
ulations. First, we take the derivative of f(x): 

_f£. 
dx 



le - xH ^|0) = B(0\e 

+ 9*Ja{x)+gtf b {x). (C13) 
Carrying a + 6 + past e~ xH yields 

~~0x = 9 * a ^ a ^ + g bh{x) + B ( uv + z i z 2)/0) + 

+Bv[zif*(x) + z 2 f;(x)} + Bv 2 (0\abe- x "\0). (C14) 
Noticing from Eq. ( p!3| ) that 

-xH,nX d f 



-xH 



a + b+\0) 



B(0\abe- xH \0) 



dx 



9af:(x)-g b f;(x), (C15) 



one obtains the differential equation on f(x) 
df 

-(1 - v 2 )— = g*f a (x) +9bh(x) + B(uv + z 1 z 2 )f(x) 
dx 

+Bv[ Zl r a (x) + z 2 fZ(x)} - v 2 lg a f:(x)+g b fZ(x)}. (C16) 



Substituting here Eq. (C12) for f a (x) and fb(x) and in- 
tegrating over x, one finally arrives at 



Tln/(/3) =A-Thx 

, {g*a+9b){ga + gt) 

2{A + B) 

, (ffa - gfeKgg - 9b) 
2(A-B) 



cosh (fie) H sinh (/3e) 

£ 



1 - 
1 - 



2T 



A + B + ecoih(f3e/2) 
2T 

A-B + ecoth(/5e/2) 



(C17) 
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